Electrical and thermoelectrical transport in Dirac fermions through a quantum dot 
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We investigate the linear conductance and the thermopower through a quantum dot between 
two electrodes of Dirac fermions with a pseudogap Anderson model using the non-crossing approx- 
imation. When the Fermi level is at the Dirac point, the conductance shows a valley, where the 
thermopower changes its sign. The conductance shows a peak when the Fermi level away from 
the Dirac point, where the impurity moment changes sharply, indicating a transition between an 
asymmetric strong coupling Kondo state and a localized moment state. The thermopower changes 
its sign in the vicinity of the peak. The magnitude of the thermopower can be more than ks/e. 
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Electron transport in graphene [1] is currently under 
active investigation. Graphene-based quantum dot (QD) 
structure is fabricated and the Coulomb blockade effect 
is observed [2-5]. In semiconductor QD systems, it is 
known that the Coulomb interaction in the QD induces 
a magnetic moment, leading to the Kondo effect [H, [?}• 
The Kondo effect in graphene has been intensively dis- 
cussed from the view of the magnetic impurity problem 
in massless Dirac fermions [8l-ll4j. A recent experiment 
shows the Kondo effect can be induced by lattice vacan- 
cies [HI . The Kondo effect is expected to be observed in 
another representative Dirac fermions, the surface states 
of the topological insulator. A single Dirac cone appears 
on the surface of three dimensional topological insulators 
[l6l-[20j. Several theoretical studies on the Kondo effect 
in this system have been developed [2ll-[24j. 

The Kondo problem in massless Dirac fermions is an 
important part of the pseudogap Kondo problem [25| . 
where the density of states of conduction electrons is 
proportional to uo r with the energy of the electrons uo 
and an index r. It is for example realized in unconven- 




2 6| | . The numerical renor- 
and perturbative scaling 



tional superconductors 
malization calculations [27|, |2£ 

analyses [29|, [30| show the system exhibits an impurity 
quantum phase transition. The massless Dirac fermion 
system corresponds to r = 1. In this case, when the 
Coulomb interaction is strong enough, three fixed points 
control the electronic states in the system, namely the 
local moment (LM) fixed point, the asymmetric strong 
coupling (ASC) or frozen impurity fixed point, and the 
valence fluctuation (VFI) fixed point located in between 
[27l . I29I . l30j . The Kondo problem in graphene has been 
studied as a tunable pseudogap Kondo problem [H, [l3| , 
where the transition can be controlled by external gate 
voltages. Recently the Kondo effect found in Ref. [15j 
has been analyzed using a pseudogap Anderson model 

& 

In addition to electric transport, thermoelectric trans- 
port in nano-structures has been actively investigated 



[321 ] . For instance, the measurement of the thermopower 
can reveal the electron-hole asymmetry in a system. The 
thermopower under the Kondo effect [33[ has been ex- 
amined in QD systems [34], [35| . In the QD system, its 
measurement clarifies the formation of the Kondo reso- 
nant state [36]. 

We investigate electric and thermoelectric transports 
through a QD between two electrodes of Dirac fermions. 
The tunneling coupling between the QD and the elec- 
trode is proportional to uj. We show how this pseudogap 
coupling controls electronic states of the system and the 
transport properties. In addition, we show how the im- 
purity quantum phase transition influences those prop- 
erties. To this end, we study the conductance, the ther- 
mopower and the impurity magnetic susceptibility based 
on a pseudogap Anderson model using the non-crossing 
approximation (NCA) (371440^. 

Model — We consider a system consisting of a QD and 
two electrodes with the Dirac fermions, as shown in the 
inset of Fig. Ufa). We treat the case of a single Dirac 
cone with the chemical potential fii(i = L, R), and focus 
on the zero bias voltage limit, //l —> Mr, setting Hi = 0. 
The position of the Dirac point from the Fermi level is 
denoted by — /io- 

The Hamiltonian of Dirac fermions in the lead z, 
is given by 



H, 



(0 



*a(fe) 



(1) 



with M{k) = ( 



), where ^ a (k) are the 
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annihilation operators of Dirac fermions, vf is the Fermi 
velocity, k = (fccos#, ks'mO) with k = \k\ and 6 stands 
for the azimuth angle of k from the x axis. The indexes 
a and b refer to (pseudo) spin indexes. For a single Dirac 
cone system, those are the spin indexes; a =t and b =|. 
For graphene, a = (A, s) and b = s) with the two 
sublattices A and which play the role of the pseudo 
spin, and the spin index s =t, I around the K and K' 
points |41j . 

The QD has a single energy level E g controlled by a 



2 



gate voltage. The Hamiltonian for the QD is given by 



H d = E z d l d s + Un t n u 
*=t4 



(2) 



where d s is the annihilation operator of an electron with 
spin s in the QD, n s = d\d s , and U is the Coulomb 
energy in the QD. The tunneling coupling between the 
QD and the lead is given by the tunneling Hamiltonian 

Ht, H T = Ei=L,R, 8 /^(^M*.(fc)+h.c.). For 
simplicity, we assume V^(fc) = V. 

It has been shown the QD with an electrode of Dirac 
fermions reduces to the pseudogap Anderson model after 
a series of linear transformations [8|, |23|, |26| . Then Ho in 
the energy space is given by 



i = L,R J 

<X=± 



duj uj c;.(w)c d (a;), 



(3) 



where c a i(uj) is the annihilation operator for an electron 
in the lead i with the pseudo spin index a at the energy 
uj. The tunneling Hamiltonian is given by 



(4) 



We have disregard the difference between the spin and 
pseudo spin indexes, both are now represented by a. 
The tunneling coupling between the dot and electrodes 
is T(uj) = a(uj + /i ) with a = 2 ^ y2 [M 0, HI . It is 
necessary to introduce an energy cut-off D, and then it 
is convenient to rewrite Y(uj) in the form: 



r( w ) 



w + fi 



D 



(5) 



with To/D = a. The pseudogap Anderson model with 
H = Hd + ffp + #0 is analyzed below. 

We consider the infinite U model. This corresponds to 
the asymmetric pseudogap Anderson model (271 [29|, [30| . 
We use the auxiliary boson technique [37], Hlj, where 
d a = b^f a , b is a boson operator and f a is a fermion 
operator; tfb + /jjg = 1 should be satisfied. Then 
we adopt the NCA 1371 kU)( to calculate the local den- 
sity of states of the QD and other quantities. To per- 
form the numerical calculations, we evaluate the Green 
functions on the range of < 10D and introduce the 
Lorentzian cut-off in the pseudogap coupling; Tj^(uj) = 
T(uj) -D 2 /(uj 2 + L> 2 ). The convergence of the NCA equa- 
tions is monitored by the unity of the spectral functions 
for f a and b . a nd the relation on the total occupancy 
in the QD [39|, which are satisfied typically with the 
accuracy of less than 0.1%. 

The conductance G and the thermopower S are given 
by G = e 2 I (T) and S = -h(T)/[eTI (T)] with 



duuj n ^^T{uj)lmA r {uj), (6) 

UUJ 
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FIG. 1. Electron transport through a quantum dot (QD) 
between two electrodes of Dirac fermions when fio = 0. (a) 
Conductance G (solid line) and occupation number in the QD 
rid (dotted line) as a function of gate voltage E g for T/To = 
1.0 x 10" 2 , and D/T = 10 (a = 0.1). Inset (a): Schematic 
picture of the QD system, (b) Local density of states Pd(w) 
at the left peak (blue), the cusp (red), and the right peak 
(green) in Fig. [TJa). 



where f(e) 
V[exp(^) 



is a 



Fermi-Dirac function, f(uj) = 
1], and A r (uj) is the retarded Green 

function for electrons in the QD [35|, |43|. The im- 
purity magnetic susceptibility x(o;) is given by xi 00 ) — 
X!^, dte < ( w+<0+ )*M (t) with M(t) = i0{t)([M(t),M(0)]), 
and M = g/J^B/^ (ft ft ~ flfi)- The static susceptibility 
X — x(0) is evaluated using NCA [38|. We set g\±B = 1 
and &b = 1 below. 

Undoped system. — Let us first consider G and the oc- 
cupation number in the QD, rid when the Fermi level is 
at the Dirac point; /io = 0. It has been shown that for 
this case, the Kondo temperature Tk = [1, [H, 125^ . 
In Fig. H{a), G and rid are plotted as a function of 
the gate voltage E g . The plot of G shows a valley; 



G ex \E a 



G m with the minimum conductance 



G m at the cusp, where E g = E*. Out of the valley, G 
decreases monotonically, consisting of double peaks. In 
the valley, rid changes gradually. This means the system 
is in the VFI regime. 

In Fig. DIb), the density of states of the QD, pd(v) = 
(-l/7r)Irm4 r (cj), is plotted. At E g = E*, pd(u) shows a 
narrow and singular peak at the Fermi level, uj = 0. The 
appearance of such a distinct peak has been discussed 
pjl, HI, 0- I n s P^e of this sharp peak of pd(w), G re- 
mains small because T(uj) suppresses the peak height. 
In Fig. El G and the thermopower S are plotted as a 
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FIG. 2. (color online) Temperature and gate voltage depen- 
dence of conductance G (a), and thermopower S (b) when 
po = 0. The gate voltage E g is scaled as (E s — E*)/T, where 
Eg is the gate voltage at the cusp of G. 
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FIG. 3. Conductance G (a), and impurity moment T\ (b) 
versus shift of the Fermi level away from the Dirac point po 
at T/To = 1.0 X 10" 4 . For a free 1/2 spin, T X = 1/4. The 
values of E g are indicated in the figures. The dotted line in 
(b) indicates T k x(Tk) = 0.07. 



function of (E g — E*)/T for several temperatures. The 
curves of G resemble one another. This means the peak 
width of the double peaks depends linearly on T. The 
curves of S are also close to each other; S cx (E g — E*)/T, 
and S changes the sign at E*. This means E* gives the 
boundary between electron- like and hole-like transport, 
This is in agreement with the fact that E* locates in the 
VFI regime. Note that the magnitude of S can be large, 
which is more than k^/e ~ 76 [/iV/K]. 

The above results of G and S are explained largely by 
the shape of pd{w). When [7 = 0, 



Pd(u) = — Im— 

7T UJ 



1 



■Eg+iT (uj/Dy 



(7) 



To adjust the QD level shift due to U , we replace E g 
with Eg — E*. When U = 0, the peak width of pd(w) 
is smaller than T, since D/Tq > 1. It is valid for the 
present model. Then pd(oj) can be approximated by a 
Dirac delta function, pd(oj) ~ S(uj — (E g — E*)), so that 
G is determined by the remaining part of the integrand in 
Eq. ©, \oj\f(uj). This factor explains G ex \E g -E* g \ and 
the width of the double peak structure. The line shape of 
S is explained as follows: h(T) cx \E g —E*\(E g —E*) near 
E g = E* g1 resulting in S cx (E g - E*)/T. The magnitude 
of S can be the order of fc#/e since f'{uo) is finite when 



\E a ~ Et 



T. 



Doped system. — Next we consider the case when the 
Fermi level is away from the Dirac point; po ^ 0. Since 



T(uj) is finite at uj = 0, the Kondo effect is active. As 
a result a clear sign of the transition between the ASC 
and LM states will appear by tuning external voltages 

In Fig. [3j G and the impurity moment T\ are plotted 
as a function of po • The values of E g are indicated in the 
figures; two of them are below E*, and two of them above 
Eg. In three of them, G shows a sharp peak. The peak 
height is about 2e 2 /h. This result is in contrast with the 
one in Fig. 0Ja). When po is negative and large, Tx is 
almost zero. This means that the impurity is screened 
by the Kondo effect. This corresponds to the ACS state. 
As po increases, Tx increases sharply, where the system 
moves to the LM state. The sharp peak of G appears 
at the boundary of the transition of those two states. 
The Kondo temperature Tk can be defined through T%; 
T K x(T K ) = 0.0701 [Hm. Here we adopt this definition 
of Tk and draw a line in Fig. [3]^b) at Tx = 0.07. Since 
Tk depends on /iq, Tk(pq) = T at the intersection of 
the line and the curves. The peak structure of G can be 
associated with the transition at T ^ Tk- Thus G(po) 
brings a clear evidence of the transition between the ASC 
and LM states. 

In Fig. |U G and S are plotted as a function of po 
when E g /Tq = —0.48 at several temperatures. As the 
temperature increases, G decreases monotonically. This 
can be understood by the suppression of the Kondo effect 
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FIG. 4. (color online) Conductance G (a), and thermopower 
S (b) versus fio for several temperatures and E g /Fo = —0.48. 
Inset (a): pd{oj) for fio — (blue) and fio = 1.2 X 10 -2 (red), 
where the dip of S appears, at T/Tq = 4.0 x 10 -3 . 



due to the temperature. In the vicinity of the peak of G, 
S changes its sign. In the ASC state, the Kondo resonant 
state locates above the Fermi level, whereas in the LM 
state, the impurity level is below the Fermi level. Thus 
the sign change of S reflects the transition between two 
states. 

We add two comments on the results of S. First, 
around fio = 0, S changes the sign for low temperatures, 
and this change disappears for high temperatures [46j . In 



this case, a dip of Pd(cj) induced by the pseudogap crosses 
the Fermi level as shown by the blue line in the inset of 
Fig. |4^a). This dip is symmetric around u = 0. The 
asymmetry of the dip becomes evident as uo is away from 
zero. At low temperatures, the crossing of the dip at the 
Fermi level results in the change the sign of S. At high 
temperatures, where the asymmetry of pd(oo) is captured, 
S remains negative. Second, around po/To = 1.2 x 10 -2 , 
a clear dip in S appears for high temperatures. This is 
the place where pd(oj) shows the critical sharp peak as in 
Fig. D^b), at uj = —p . 

We have shown that S can be more than /ce/e, which 
is dissimilar to the conventional Kondo system (38^ . 
The figure of merit, ZT = S 2 GT/n, where k is the 
thermal conductivity, describes the efficiency of thermo- 
electric materials |32j, and it can be calculated using 
the Landauer formula approach, where n = [h(T) — 
h(T) 2 /I (T)]/T [H El- Within the present model, 
ZT > 1, for wide range of parameter values of T, E g , 
and /io, and sometimes ZT becomes more than 10. Thus 
the pseudogap Kondo system is interesting not only from 
the point of view of the impurity quantum phase transi- 
tion, but also for its thermoelectric properties. 

In summary, we have investigated the conductance and 
the thermopower through a QD between Dirac fermions. 
The pseudogap tunneling coupling between the QD and 
electrodes suppresses the conductance for undoped elec- 
trodes. For doped electrodes, the conductance shows a 
peak structure which signifies the transition between the 
ASC and LM states. The thermopower can be more than 
k^/e for both doped and undoped cases. 
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